General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



1 







) 


'I 

I ■. 

-1 


r 

I. 


I 

EFFECTS OF STREM- ASSOCIATED FLUCTUATIONS 
UPON THE RADIAL VARIATION OF 
average solar-wind PARAMETERS 



' i • 

! '! ' 


O 


I 

.■f 









iS .6 




by 

Bruce E. Goldstein* 
J. R. Jokipiit 


*Jet Propulsion Laboratory, Pasadeng., California 91103 
Temporarily at National Aeronautics and Space Administration 
Code SL, Washington, DC 20546 until August, 1977 

tDepartment of Planetary Science? , Department of Astronomy 
Unviersity of Arizona, Tucson, Arizona 85721 


ABSTRACT 


The effects of nonlinear fluctuations due to solar-wind streams 

upon radial gradients of average solar-wind parameters are computed, 

using a numerical MHD model for both spherically symmetric-time 

dependent and corotating equatorial -flow approximations. We find 

significant effects of correlations between fluctuations upon the 

gradients of azimuthal magnetic field, radial velocity, density and 

azimuthal velocity. The averages of and r^n both 

decrease with increasing radial distance (O to 30 km/sec, 5 to 104 

by 6 AU) while mass flux is conserved. The cross-variance between 

V and r^n increases with radial distance from negative values to 
r 

positive values at distances greater than 3.5 AU. The average azi- 
muthal velocity decreases , and at distances greater than 4 AU may be 
opposite to the solar rotation ; <V^> of from 0.2 to -1.7 km/sec at 
6 AU is predicted. Between ^00 to 900 solar radii stream inceractions 
have transferred the major portion of the angular-momentum flux to the 
magnetic field; at even greater distances the plasma again carries 
the bulk of the angular-momentum flux. The average azimuthal com- 
ponent of the magnetic field may decrease as much as 10^ faster than 
the Archimedean spiral out to 6 AU due to stream interactions, but 
this result is dependent upon inner boundary conditions. 
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1. INTRODUCTION 

Recent observations carried out over large ranges of helio- 
centric distances with deep space probes have determined the radial 
gradients of various solar-wind parameters (velocity, magnetic field, 
etc.). Theoretical analyses of the radial variations have conven- 
tionally been carried out by attempting to construct steady-state 
fluid dynamical models which reproduce the time averaged parameters. 
However, this approach neglects possible effects of correlations 
between large fluctuations in the parameters; the observed, large, 
non-linear fluctuations can be correlated to change the average 

properties of the solar wind (Siscoe, 1970; Jokipii, 1975, 197^) in 
a way that cannot be reproduced in a non-fluctuating model. For 
example, Schubert and Coleman (1968) have suggested that MHD fast 
waves may carry a substantial portion of the solar wind angular 
momentum fl ux, Pizzo et al . (1973) and Neugebauer (1976) have shown 
that fluctuations will alter the expected relation between temperature 
and velocity in the solar wind, and Jokipii (1975) has suggested fluc- 
tuations as a possible explanation for the faster than 1/r decrease 
in the azimuthal interplanetary magnetic field (but Parker and 
Jokipii (1976) show that for Pioneer 10 data Inclusion of the varia- 
tion of the average solar wind velocity considerably improves 
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agreement with the Parker spiral field). 

In the theoretical analysis (Siscoe, 1970; Jokipii, 1975, 1976), 
the relevant fluctuating quantities which were shown to affect the 
mean were the correlations between fluctuations, and their radial 
variations. The magnitudes of these correlations were left either 
to direct measurements or to theoretical calculations based on 
solai — wind models. It is the purpose of this paper to present the 
results of detailed numerical calculations intended to reveal the 
existence and possible magnitude of the effects of fluctuations on 
the averages of selected solar-wind parameters. Since variations and 
correlations associated with stream structure at periods greater 
than one day are observed to dominate the fluctuations in the solar 
wind (Goldstein and Siscoe, 1972), the calculations presented here 
concentrate on the effects of stream-related phenomena on the radial 
dependence of average solar^-wind parameters. 

The organization of this paper is as follows: In the next 

section the model and computational techniques are described in 
detail. Section 3 discusses effects of stream structure upon the 
radial gradient of the density-velocity correlation in the solar 
wind, and the consequences for deceleration of the solar wind and 
estimates of solar vrind heating. Section 4 continues this dis- 
cussion in the context of angular-momentum transport and the average 
azimuthal velocity; it is demonstrated that at distances of h AU and 
beyond the average azimuthal velocity may be in the direction 


oppos i te to the solar rotation. In section 5 the variation of 
with radius is Investigated; for certain sets of inner boundary 
conditions (not the most realistic) the effect of stream variations 

is to distort the average of B, from the Archimedean spiral. In 

<P 

section 6 we summarize our results and discuss the observational 


consequences . 
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2. EQUATIONS, BOUNDARY CONDITIONS, AND NUMERICAL SOLUTIONS 

Our starting point is the equations of motion for an MHO fluid 
written in conservation-law form (Rossi and Olbert, 1970). Ignoring 
displacement currents and relativistic effects, these equations can 
be written as follows: 


l2- 

8t 


+ V-dV = 0 


( 1 ) 


I- ?V = -vCpVV + P - 1/Uo(BB - iB^Sik)] + pg 

O t ^ 


( 2 ) 


= 7x(Vx^ 
3t 


(3) 


V-B = 0 




1^ + e + BV2y^) = 


-7-[(ipV^ + e)V + P-V + q + 1 /Uq Bx(VxB)] + pg-V 


(5) 


In these equations p is mass density, V is solar wind velocivy, B is 
magnetic field, e is internal energy density of the plasma, P is 
the plasma pressure tensor, q is the heat flux, and g is the solar 
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gravitatlonal acceleration. For the purposes of solving these equa- 
tions we adopt a coordinate system in which r is radial from the sun, 


9 is the longitudinal coordinate, and 0 is colati tude^and a model for 
plasma flow in the solar-wind equatorial plane in which Bg = = 0 

and — = 0. The limitations of this approximation (which ignores 

39 3V 3B 

meri dional flows, i.e. and terms) are discussed in Suess 

and Nerney (1973). The primary effect is that flow divergence from 

the equatorial plane due to the r sin 0 dependence of is ignored, 

Nerney and Suess (1975) snd Winge and Coleman (197^) show that this 

procedure does not cause decreases of from more than 10 to 204 in 

B at 5 AU due to flow divergence; somewhat larger effects may occur 
6 

for V . Rhodes and Smith (1975) have presented observational evi- 
dence for an increase in solar wind velocity with latitude; gg terms 
will also effect, radial gradients in the equatorial plane. The 
effects we discuss will also occur in models that include latitudinal 
variations; such effects can be regarded as contributing in a some- 
what independent fashion. Finally, we note that at high lati- 
tudes v/ill be less than in the equatorial plane due to the sin 6 
dependence of rotational velocity; Suess et al. (1975) discuss 
theoretical models of latitudinal variations of solar wind streams. 

In the equatorial-flow approximation we have 





( 6 ) 
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( 7 ) 


3 r 


- V B ) = 

(I) r r (p 



( 8 ) 


Ip r^LpV/ + i(B ^2 - B^ 2 )/y^ P] = 

1^ r(B^B/„„ - pV^V^) - Ip r^pV^ 

+ r( 2 P + p(V^^ - G/r) + (9) 

|pr3(pV^V^ - B^B^/y„) . 

|_ ,a(.p . . |_ (,0) 

Ip r^[V_.(ip(V^=^ + V^ 2 ) + P + £) + 1 /u„ - V^B^)] = -pGVj. 

^ ^ C-P%(ip(V,^ H- ^ p) ^ ,Bp(B^V^ - V^B^)/y„] 

- Ip r2[Up(v/ + V^2) + + (B^2 + b^2)/2u„] (11) 


( 12 ) 
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In eq. 11 we have assumed that heat flux Is zero; and we assume 
that pressure is a scalar, P. The effects of thermal conduction are 
approximated by a polytrope relation between the internal energy 
density and the pressure of the plasma. Thus, 


P = P + P. 
e I 






1 


P . 

e, I 



(13) 

(1M 

(15) 

( 16 ) 

(17) 


In these equations a is a constant mul tipi ier for the polytrope law 
(eq. 1^) that is conserved along stream lines (eq. 17). The sub- 
script e or i refers to the pressure, internal energy, and polytrope 
index of the electron or Jon gas. The proton temperature is deter- 
mined by conservation of energy rather than a polytrope law. This 
scheme was introduced for the purpose of another study; the results 
to be described do not depend critically upon the particular tem- 
perature law assumed. 

In this study we investigate both spherical ly-symmetric, 
time-dependent solutions of the solar-wind equations, and longitudinally 
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varying solutions that are time independent in a framework corotating 
with the sun. For the spherically-symmetric cases, we set “ in 
eqs. 6-12; whereas for the corotating case the time derivative in 
inertial space is replaced by where P. is the solar angular 

velocity. Inner boundary conditions are then specified at r - 0.2 AU 
as a function of t or <{.. and the MHD equations are integrated numer- 
ically outward to 6.0 AU. The longitudinal boundary conditions 
are periodic in <|); the time-dependent boundary condition is chosen 
to be periodic so that an analogous computation is possible. 

For the spherically-symmetric model, eq. 7 is simply the state- 
ment that r^B^ is conserved, but eq. 8 must be integrated numerically 
For the corotating model eq. 7 must be integrated numerically, but 
it is possible to avoid integrating eq. 8 by choosing an initial 
condition in which fluid velocity and magnetic field are aligned. 
Thus, eq. 12 requires that 


I* [Vr “ 


(18) 


in the corotating case; we actually determine at the inner boundary 
by the relation 





( 19 ) 
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after other variables have been specified. As the electric field 
is curl-free, eq. 19 is satisfied everywhere (not just on the Inner 
boundary), and is a conserved quantity. 

To facilitate the numerical calculation, we adopt a change of 
variables; 


Vq 


/ O*. J. 

or A = -fit + — 
o 


(20a, b) 


where V is the average solar-wind velocity. The radial variable is 
o 

not altered and the transformation replaces the original longitudinal 
variable with a longitudinal variable that remains constant along 
the average spiral direction. The change of variables introduces 
additional derivatives with respect to ' on the right hand side of 
eqs.6-12 but does not affect the conservation law form. This change 
of variables allows a larger step size in Ar. 

Equations 6-12 can be written in vector form; 


8A _ 8G , 

IF " 


( 21 ) 


where A is a conserved quantity of the flow for those equations for 
which the I component Is zero. The numerical calculations are con- 
ducted using a predictor-corrector technique. Values of A at r + Ar 
are calculated from the right hand side of eq. 21 evaluated at r. 
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A at r, A at r - Ar, and a finite difference operator. The differ- 
ence operator is an adjustable combination of a second-order 
accurate difference operator and a small contribution from a first- 
order accurate difference operator chosen to provide numerical dissi- 
pation by averaging ^ in the ^ direction. The plasma parameters at 
r + Ar are calculated from A and the auxiliary equations by a 
rapidly convergent iterative technique. The value at r + Ar is then 
corrected by use of eq. 21 with derivatives and source terms evaluated 
at r + Ar. 

The calculation used a variable step size, Ar, restricted 
to satisfy the Courant condition (Fox, 1362). In addition, a 
limit Ar < 0.02r was imposed; in practice this condition determined 
the step size. In the variable 128 steps were used and the period 
of repetition v/as chosen to be 1/3 of a solar rotation for both 
corotating and spnerically symmetric solutions. The stability and 
accuracy of the solutions were investigated using a variety of 
different step sizes In the r and ^ variables and by altering the 
amount of numerical dissipation. For the calculations to be 
described there are no significant errors in the average properties 
of the flow parameters; further decreases In step size only improve 
resolution in the vicinity of shocks that form in the flow. 

To determine the effects of solar-v^Ind stream structure upon 
the average properties of the solar wind, we have performed calcula 
tions using a variety of inner boundary conditions (both realistic 
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and unrealistic) to aid in understanding what effects are important 
in causing the radial gradients that occur in the solution. The 
results we discuss in this paper are based on the following set of 
inner boundary conditions Imposed at 0.2 AU: 

= ^50 + 150 (cos 3^), (km/sec) (22a) 

M ^ - 450 + 200 (cos 3 X) (km/sec) (22b) 

V, = 40 (km/sec) (23) 

9 

n = 175 , n - 175 (l - i cos 3/), (n/cm"^) (24a, b) 

n = 175 (1 - .65 cos 3X) (24c) 

T = T. = 3.3x105 (®K) ( 25 ) 

e f 

= 100, = 100 (1 - i cos 3^), (y) (26a, b) 

= 125 (26c) 
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Throughout the paper we will discuss primarily calculations 
using eqs. 22a, 23, 24b, 25, 26a as inner boundary conditions (Case l) 
to determine what happens in the absence of density-velocity corre- 
lations we substitute 24a for 24b (Case 2); to consider large 
an t i corre 1 a t i ons we use 22b, 23, 24c, 25, 26c (Case 3), and in the 
section on azimuthal magnetic fields we also use 26 b (rather than 
26a), 22a, 23, 24b, and 25 (Case 4). These choices are summarized 
in Table 1, where the entries refer to equation numbers. 

A polytrope Index of 1.2 Is chosen for the electron gas and 
5/3 for the ion gas (the proton polytrope index affects only the 
relation between internal energy and pressure because the proton 
temperature is determined by the energy conservation equation). 

Note that our procedure is different from that of some other 
workers (e.g., MHD model of Steinolfson et a I . (1975)) who specify 
conditions for ail r at t = 0 and then Integrate forward In time. 

We specify a cyclic boundary condition for all time at r = 0.2 AU 
and then integrate outward in radial distance. At larger distances 
the evolution of our solutions Is similar to that shown in Fig. 2 
of Hundhausen and Gosl Ing (1976) ; the velocity gradients In regions 
between overtaking streams steepen, forward and reverse shocks form 
on either side of the compression region, etc. The averages which 
we report in the following sections are obtained by integrating the 
solar wind parameters for such solutions over a periodicity cycle, 
at a fixed radial distance. 
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3. DENSITY-VELOCITY CORRELATIONS AND RADIAL GRADIENTS OF PLASMA PROPERTIES 

A salient feature of the solar wind which was first noted by 
Neugebauer and Snyder (1966) , and by numerous workers since that 
time, is the existence of high-speed streams in the solar wind and 
an anticorrelation between radial velocity and density. As high- 
speed streams overtake low speed plasma, a compression region 
develops along the interface between high and low-speed streams. 

In the case of a source region that rotates with the solar corona, 
this interface will lie along the average spiral direction of the 
interplanetary magnetic field. (In the case of a spherically 
symmetric time-dependent model compression will occur in a spherical 

shel 1 . ) 

To consider how such an interaction might evolve with distance, 
we have integrated the MHD equations of motion for periodic inner 
boundary conditions at 0.2 AU (Case 1, previous section) with a density 
and velocity anticorrelation; for the purposes of comparison we have 
also investigated conditions (Case 2) with high speed streafij'S but 
with an initially uniform density distribution, and conditions 
(Case 3) with a larger than average stream structure and associated 
density-velocity anticorrelation. The numerical calculations were 
done for both spherically symmetric-time dependent models and 
time-independent, corotating, azimuthally varying models. Tne 
value of radial velocity, averaged over a periodicity cycle (time 
or longitude), is shown in Fig. 1 as a function of radial distance. 
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It can be seen that the effects of velocity-density correlations 
are much greater than the differences between symmetric and 
corotating models. The predicted decrease in velocity with radial 
distance is due simply to conservation of momentum; if a high speed, 
low density object collides with a low speed, high density object, 
the low density object will be more greatly decelerated than the high 
density object is accelerated. This is a non-negl igible effect; 
conversion of internal energy and heat conduction energy into 
kinetic energy at large radial distances v/ould lead to an expected 
increase of radial velocity of about 25 km/sec between 1 and 6 AU, 
but actual decreases of 7 to 25 km/sec are calculated. The differ- 
ence of about 30 to 50 km/sec is due to the non-linear effects of 
the solar v/ind density-velocity correlation upon the equations for 
average solai wind flow. 

The consequences of this decrease in averaged with increasing 
radial distance for the radial gradient of solar-wind density can 
be determined by averaging eq. 6 (flux conservation) over a period- 
icity cycle (longitude or time). We then obtain 

C = r^<nV^> = + <5n6V^>) , (2?) 

Vi'here c denotes difference of a variable from its average value, 
brackets denote an average value over longitude (or time), and C is a 
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constant. We have plotted r^<n><V^> and r^<6n5V^> in Fig. 2 for a 

spherically symmetric integration of the Case 1 boundary conditions. 

The ant icorrelation between V and n decreases with increasing 

r 

radial distance, and in fact density and velocity are positively 
correlated at distances greater than 500 solar radii. This rather 
surprising result arises from a rebound phenomenon; If a high 
speed, low density stream collides elastically with a low speed, 
high density stream, the velocity of each stream will simply be 
reversed in the framework moving with the center of mass. Thus, 
density and velocity will be positively correlated in the compression 
region between forward and reverse corotating shocks that form in 
the solar wind; at large radial distances these regions expand so 

that most, and eventually all, of the solar wind has been shocked. 

The consequences of this phenomena for radial gradients of density 
are significant; the increase in the r^<6n6V^> term implies a corres- 
ponding decrease in r^<n> with radial distance (6% from 1 AU to 
6 AU) despite the fact that is approximately constant. Thus, a 
steeper than 1/r^ fall off of <n> need not be interpreted as evi- 
dence for a meridionally diverging flow. 

To further elucidate the importance of fluctuations we have 
plotted in Fig. 3 the terms that contribute to radial gradients of 
radial momentum flux for the same calculation (spherically symmetric 
Case 1). (Radial momentum flux can not be written in strict conser- 
vation law form in a spherical geometry due to geometric source termn.) 
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Contributions to gradients of <nxV^>2 at distances of llO 
to 160 solar radii from both pressure and nonlinear fluctuations are 
significant; at distances greater than 160 solar radii the fluctua- 
tion term dominates. In fact, a decrease of about 15^ in the 
quantity r^<n><V^>^ occurs beyond 160 solar radii, while the actual 
momentum flux carried by the bulk flow undergoes a very slight 
increase with distance. Thus, a decrease of solar-wind velocity with 
increasing heliocentric distance must be interpreted with caution. 
Stream-stream interactions and the resultant correlation between 5n 
and can explain a velocity decrease without necessarily requiring 

energy or momentum sinks such as pickup of interstellar neutrals or 

acceleration of energetic particles. 

We have also computed the radial variation of kinetic energy 
flux (not shown, but dominated by fluctuations similar to Fig. 3); 
a minimum value of this quantity v^<nM Vl> occurs at a model depen- 
dent distance of ^30 solar radii. At this distance the kinetic 
energy flux initially present in the form of stream kinetic energy 
(fluctuations) has been converted into internal energy of the flow 
(gas and magnetic). At greater radial distances the flow expands, 
the energy Is reconverted into a kinetic flux, and the proportion of 
energy flux in the bulk flow increases. This interaction is reflected 
In the radial variation of gas and magnetic pressure (Fig. 3) where 
local maxima of values integrated over a spherical shell (scaled by 
r^) occur at distances of solar radii. 

An upper limit for the effect of stream-stream interactions 
upon the heating of protons can also be estimated. Ion-temperature 
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gradients might be due to either of two separate effects: adiabatic 

temperature variations due to compression and expansion of the 
plasma, and dissipative heating due to corotating shocks. The 
extent of dissipative heating will depend upon the relative propor- 
tions of energy at shocks that is taken up by the electron component 
of the gas. The effects of electron heat conduction (see Section 2) 
are approximated in the case of a poly trope law by an increased 
internal energy (eq. 16) for a given temperature. In addition, the 
macroscopic flow characteristics such as distance of shock formation 
and amount of compression should also depend upon radial gradients 
of electron temperature (and hence polytrope index). Given that 
some assumptions concerning an electron heat law are made, the 
proton temperature can then be computed using conservation of energy 
flux (eqs .11,16). 

Proton temperature profiles for electron polytrope indexes of 
1.2, 1.5, and 5/3 for Case 1. inner boundary conditions with spherical 
symmetry are shown in Fig. Surprisingly, there is very little 
difference between the temperatures predicted for electron polytrope 
indexes from 1.2 to 5/3. We have also plotted electron-temperature 
gradient for the case = 5/3- As our model puts all dissipative 
heating in protons, the temperature profile shown represents solely 
the effects of adiabatic heating by compression. It is apparent that 
the effects of plasma-density variations upon average values of adia- 
batically determined temperatures are negligible; the gradient follows 
closely the power law that would be expected for a model with- 

out temporal (or longitudinal) variations. 
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A1 though our treatment of the electron heat conduction problem 
Is far from ideal, the wide range of electron poly- 

trope indices used should bracket reasonable models for the behavior 
of the electron gas. Our calculations show a significant minimum in 
the solar wind proton temperature at 300 solar radii; In the actual 
solar wind this minimum will be more spread out due to the variable length 
scales of actual solar wind streams and may not occur at all. A broad 
plateau v/i th temperatures on the order of 8x10^°K at distances of 
from 3 to 6 AU is a feature that should be readily apparent in solar 
wind data. However, this is not guaranteed. For the case of elec- 
tron polytropic index y = 1.2, the internal energy (longitudinal 
average of p ) wi 1 1 not be greatly affected by longitudinal density 
gradients. For y^ = 5/3, a substantial percentage increase in the 
internal energy of the electron gas may occur. However, electron 
temperatures at 1 AU and beyond are unrealistically low due to the 
inner boundary condition assumed (it is not possible to model accu- 
rately temperature with a polytrope law over a wide range of radial 
distances). Thus, more energy may be absorbed by adiabatic com- 
pression of an electron gas than indicated by the calculations. In 
addition, dissipative heating of the electron gas as well as the 
ion gas is to be expected. The ion temperature profiles calculated, 
although they might be realistic, are best interpreted as plausible 
upper limits on heating due to stream-stream interactions. A broad 
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plateau with temperatures on the order of 8 xlo’^°K at distances of 
from 3-6 AU is suggested, and may be apparent in the solar wind. 

An important point to note is that the large decrease in the standard 
deviation of with increasing radial distance does not imply 
counter balancing increases in thermal speed. As discussed earlier, 
transfer between terms in the energy flux equation is much less than 
might be expected from gradients of products of averaged density 
and velocity. 
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k. ANGULAR'MOMENTUM TRANSPORT AND AZIMUTHAL VELOCITY 

Measurements of the solar~wind azimuthal velocity and angular- 
momentum flux are reported by Lazarus and Goldstein (1971) 

from Mariner 5 plasma and magnetic data; Egidi and Signorinl (197^) 
present measurements of azimuthal velocity obtained from HEOS 1 and 
Explorer 33 spacecraft; Hundhausen et al . (1970) report Vela 3A and 
3B measurements; and EgidK^^. (1969) give Explorer l8 results. The 
concensus of these observations, despite some disagreement, Is an 
average V^ at 1 AU on the order of 7 to 8 km/sec. These measure- 
ments are, however, difficult to obtain, and other effects such as 
variations on periods of months to years (Egid; and Signorini, 197^) 
may also be occurring. Tangential velocities as large as 8 kra/sec 
are larger than vvould be expected on theoretical grounds (Weber and 
Davis, 1967 ), but theoretical uncertainties (Weber, 1970, 1972) and 
non-uniform solar-wind expansion (Priest and Pneuman, 197^) are also 
significant. Further discussion of observational and theoretical 
problems are provided by Hundhausen (1972) and in Solar Wind (edited 

by Sonett aj_. , 1 972 ^ pages 26J--275). 

For the purpose of estimating radial gradients beyond 1 AU, an 

inner boundary condition of V^ = ^0 km/sec Is chosen at 0.2 AU; this 

results in an average V, of 8 km/sec at 1 AU. Using a radial mag- 

netic field at o)2 AU that results i n a B . of at 1 AU, 

(j) 

approximately 2/3 of the angular momentum flux at 1 AU is carried by 
the plasma and 1/3 by the magnetic field; these values are in 
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approximate agreement with the 3 to 1 ratio observed by Mariner 5 
(Lazarus and Goldstein, 1971). As conditions match observations at 
1 AU, the predictions for angular- momentum transport and azimuthal 
velocity at larger radial distances can be directly determined and 
compared to observations. 

Unlike other plasma parameters, radial gradients of azimuthal 
velocity and angular momentum flux depend sensitively upon whether 
a spherically symmetric or corotating model is used. The physical 
basis for this difference is the orientation of the high pressure 
compression region between low and high speed streams. In the 
spherically symmetric model this compression region is a spherical 
shell, v/hile in the corotating situation the high pressure region is 
aligned along a spiral interface. ThuSj in the spherical ly-symmetri c 
case pressure gradients act only in the radial direction, while in 
the corotating case pressure forces also deflect low-speed streams 
in the direction of solar rotation and high-speed streams in the 
opposite direction. The effects of this difference can be seen in 
radial gradients of radial velocity (Fig. l). Recalling that the 
average radial velocity decreases because low-density, high-speed, 
regions are more accelerated than high-density solar wind, the 
reduced radial pressure gradients of the corotating case result in 
smaller radial accelerations and less decrease in the average radial 
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velocity v/ith increasing distance. Another effect of the spiral 
alignment is that high-speed streams are deflected against the 
direction of solar rotation and low-speed regions are deflected in 
the direction of solar rotation (Si scoe et al . , 1969). A further 
point that has not been noted in this regard is that the low mass flux 
(high speed) regions are more greatly deflected than the high mass 
flux (low speed) regions in the azimuthal direction with a con- 
sequent decrease in average a2imuthal velocity with increasing 
radial distance below the value that would be expected for a solar 
wind without fluctuations. In fact, we shall demonstrate that at 
large radial distances (900 solar radii or more) the average solar- 
wind azimuthal velocity can be in a direction opposite to the solar 
rotation . 

We investigate the effects of density-radial velocity 

correlations upon azimuthal velocity gradients by integrating models 
with and without such correlations (Case 1 and Case 2) for spheri- 
cally symmetric and corotating models (see Fig. 5). The calculations 

are started at 0.2 AU with rV^ of 8 AU km/sec. For models with no 

<P 

magnetic field or no non-linear fluctuations this quantity would be 
constant with radial distance if pressure goes to zero sufficiently 
fast as r Weber and Davis (196?) have shown that the 

effect of including magnetic fields in a spherical ly-symmetric 
model results in a transfer of angular momentum from the magnetic 
field to the plasma with increasing radial distances to yield 
slight increases in r<V.> wi th d istance. In the calculations 


vie show there is a slight decrease in r<V,> between ^0 and 70 solar 

radii; v;e attribute this to the fact that our solutions do not match 

properly through the Alfvenic and sonic transition points from a 

flow that originates in the corona. The decrease near the inner 

boundary therefore corresponds to damping of inward propagating 

Alfven waves closer to the sun. From 70 to A25 solar radii there is 

an increase of about 20% in r<V,> for the solutions shown; this 

<P 

modest increase does not occur in our models unless we include both 
the effects of magnetic fields and solar wind stream structure. The 
increase is not due to a transfer of angular momentum from the field 
to the plasma^as our later discussion will show. A similar small 
increase in the average value of i*<V^> between the corona and 1 AU 
was demonstrated in an MHD model that included stream'"Stream 
interactions (Urch, 1972). 

At larger radial distances (where reexpan^-ion of the compressed 
plasma dominates) a much larger decrease in the average value of 
r<V > is predicted. This decrease is due to the deflection of low 
density, high speed streams in a direction counter to solar rotation. 
High density plasma ahead of the spiral compression region is 
accelerated in the corotating direction, but as the density is higher 
while the pressure forces are equal, the acceleration is less. 
Consequently, a net decrease in the value of t<V^> with distance 
occurs without necessarily altering the proportion of angular momen- 
tum carried by the plasma. 
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Surprisingly, a similar phenomenon can occur in spherically- 
symmetric flows. In Fig. 5 , the spherical ly-symmetric calculation 
with a density-velocity anticorrelation (Case l) also shows a decrease 

In r<V.> with radial distance, although the gradient is only about 

<? 

half as steep as in the corotating case. The decrease may be 
attributed to magnetic forces; the bending of field lines in the 
vicinity of a compression region will tend to accelerate fast plasma 
counter to the solar rotation direction and slow plasma in the solar 
rotation direction. As in the corotating case, the density differ- 
ence between the two streams results in greater acceleration of the 
high-speed stream with a consequent decrease in the average value 

of rV. with radial distance. 

9 

An interesting question suggested by these results is whether 
or not the stream-stream interactions also cause a transfer of 
angular momentum flux between the magnetic field and the plasma. 

This question is answered by the results shown in Figs. 6 and 7. 

The average plasma angular- momentum flux, normalized with respect to 
an expected r”^ dependence, (r/AU) ^<nV^V^> (protons/m sec^) is plotted 
as a heavy solid line; brackets denote averages over longitude or 
time. The angular-momentum flux carried by the magnetic field, 

-Cr/AU)3<B^B^ >/y^mp, IS plotted as a thin solid line; the angular 

momentum fluxes due to fluctuations, -(r/AU)^(<nV V,> - <n><V ><V^>), 

and (r/AU)3(-<B B > + <B ><B >)/u m , are shown as heavy and thin 
r ({) r (J) o p 

dashed lines. For the case of the spherical ly-symmetric model this 
latter quantity is identically zero as cannot vary with time. 
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Interpretation of the spherical ly-symmetric solution (Fig. 6 ) 
is straightforward. No significant transfer of angular momentum 
between the piasma and magnetic field occurs at large radial dis- 
tances; this is in agreement with the conclusions of Weber and 
Davis 01967). The differences between the anguiar momentum flux 
carried by the plasma and the estimate that might be obtained by 
multiplying appropriate average values of density, radial velocity, 
and azimuthal velocity together are significant. A comparison of 
the average value of rV^ with radial distance (Fig. 5 ) to the 
difference between averages of products and product of averages 
(dashed heavy line of Fig. 6 ) demonstrates that the difference can 
be attributed di rectly to correlations between rV^ and r^V^ that 
develop as a result of compressions in the solar wind. 

The angul ar— momentum flux for the corotating calculations is 
shown in Fig. 7 . This case is significantly different in that a 
iarge transfer of angular momentum from the plasma to the magnetic 
field occurs as radial distance increases from 200 to 6OO solar 
radii. This can be readily explained by the compression of the 
magnetic field. The thin dashed line indicating momentum flux 
carried by magnetic fluctuations exactly follows the increase in 
total angular momentum flux carried by the plasma. This is in 
agreement with the idea that transfer of angular-momentum flux to 
the average magnetic field must be quite small; elsewhere in this 
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study we show thats.Bp>and<B^do not vary greatly from r”^ and r”^ 
power law behavior as a result of stream-stream interactions. 
However, compression of the magnetic field into the region between 
corotating streams results in a non-linear Increase in the average 
value of the product although the average values of B^ and B^ 

are not increased. The difference between the corotating calcula- 
tion and the spherical ly-symmetric calculation occurs because is 
increased by compression in spherically symmetric models. At dis- 
tances from ^00 to 600 solar radii, kinetic compressive forces are 
approximately balanced by expansive forces due to gas pressure and 
the magnetic field. At radial distances greater than 600 solar radi 
the angular-momentum flux carried by the magnetic field begins to 
decrease. At these distances, the velocity differences across 
interaction regions have decayed considerably (Collard and Wolfe, 
197^; Smith and Wolfe, 1976; Gos 1 i ng ^ al_. , 1976). Thus, the 
rebounding plasma expands away from the compression region and the 
maximum values of <( r^B^> and<rB^>thereafter decrease with 
increasing radial distance. 


not 


5. VARIATION OF AVERAGE AZIMUTHAL MAGNETIC FIELD WITH RADIUS 

The radial variation of the azimuthal interplanetary magnetic 
field B is a subject of considerable observational and theoretical 
interest. A number of observations have been reported of the radial 
variation of the average of (Burlaga and Ness, 1968; Coleman and 
Rosenburg, 1968; Rosenburg, 1970; Coleman et a1 . , 1969; Rosenburg 
and Coleman, 1973; Smith, 197^; Villanti and Mariani, 1975; 

Behannon, 1975). The observations have generally indicated that the 
average of (averaged over time periods of the order of 27 days) 
tends to fall off more steeply than the 1/r dependence expected for 
a constant solar wind. Jokipii (1975, 1976) pointed out that short 
period fluctuations in the wind velocity could change the radial 
dependence of averaged, solar-wind parameters. In this section we 
examine the numerical solutions described in section 2 and show that 
for certain sets of boundary conditions, the effect of fluctuations on <B 
can be significant. 

Before discussing the numerical calculations, we discuss the 
expected radial variation of the average B^ in terms of averages and 
correlations between "fluctuations. The evolution of the (j) component 
magnetic field is determined by eq. 8. 

We suppose that the wind velocity is the superposition of a 
mean value plus a fluctuating part 
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Slmilarly, we set B = <B> + SB. Taking the average over <f> of eq. 8, 
the derivative with respect to (j) drops out, depends only on r, 

and we can define the variable by 

A = <V ><B >-<VB>=“- <6V 5B,> + <5V,5B > (29) 

I r ({) (j> r r r 4) <1> r 

eq. 29 is a slight generalization of eq. 5b of Jokipii (1975) J C is 
an arbitrary constant. As <V^>, <B^>, and Be determined 

from other equations, can be found from Note that in the 

absence of fluctuations, and since 's small compared with 

one expects that if <V^> is constant, <B^> “ 1/r, as in the Archi- 
medean spiral. it is a general consequence of V-B = 0 that 
r^<B^> is a constant, even in the presence of arbitrarily large 
fluctuations. Thus eq. (29) in the absence of fluctuations corres- 
ponds to the conventional spiral magnetic field. Deviations of A^^ 
from a 1/r dependence imply deviations from the spiral. 

The computer-generated solutions were examined to determine the 
relative importance of the various terms appearing in eq. 29. The 
average values of V and B were constructed at each radius by averaging 
over solar longitude. The correlations of the fluctuations were con- 
structed by computing the difference between the value at each 
point and the average, and then computing the longitudinal average 
of the relevant products of the differences. As the numerical 
integration i s done in conservation law form, eq. 29 is exactly 
satisfied by the solutions obtained and an exact r”^ dependence for 

B also resul ts . 
r 
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We found In our corotating model that with Inner boundary condi- 
tions including a density-radial velocity anticorrelation and constant 
magnetic field (Case l), the left hand side of eq. 29 could be 
accurately represented by a 1/r dependence; so, the effect of 
fluctuations was not great. However, without the density-velocity 
ant i correl at ion (Case 2), or with an inner-boundary correlation 
between magnetic field and velocity (Case 4), the effects of the 
fluctuations were significant. Shown in Fig. 8 is the quantity 
A = <V ><B,> ^ <V.xB > plotted as a function of r for Case 4. For 
comparison a 1/r curve is also shown. Clearly departs signi- 
ficantly from a 1/r dependence. A least squares fit to between 
1 and 6 AU gives A, (r) = r”^*°^. As the term <6B 5V > is much 
larger than <6B^6V^>, fluctuations act upon B^ to distort the spiral 
interplanetary magnetic field. To demonstrate this we have evalu- 
ated the right side of eq. 29, including the correlation term. The 
quantity A^ + <5V^5B^> = A^ should vary as 1/r, since <6V6B^> is 
small. The quantity A^ has been evaluated in this case and follows 
the expected 1/r dependence quite well. A least squares fit to A^ 
gives a power law dependence of 

We conclude that with inner boundary conditions corresponding 
to Cases 2 and k the average azimuthal magnetic field can deviate 
significantly from the Archimedean spiral based on the average 
radial wind velocity. ' 
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discussion AND SUMMARY 

There has already been a good deal of discussion of the role of 
St ream- stream Interactions in reducing the amplitude of velocity varia 
tions in the solar wind at large radial distances (Collard and Wolfe, 
197 ^; Smith and Wolfe, 1976; Hundhausen and Gosling, 1976; Gosling 
et al_. , 1976 ). In addition, our calculation predicts a modest 

decrease of the solar wind speed at large radial distances. Collard 
e^ al_. ( 1975 ) have reported a tendency for the solar wind velocity 

between Pioneer 11 and Pioneer 10 to decrease when velocities at 
Pioneer 11 are above 450 km/sec; In addition, a decrease between 
Pioneer 11 and Pioneer 10 of the average velocity of about 30 km/sec 
occurs (Wolfe, personal communication). Our model calculations pre- 
dict a decrease from 1 to 6 AU on the order of 7 to 25 km/sec, as 
opposed to the 25 km/sec increase that would be expected from con- 
version of thermal and magnetic energy into kinetic energy in models 
without fluctuations. The magnitude of the deceleration that we 
predict depends obviously upon the magnitude of initial density- 
velocity anticorrelation assumed at 0.2 AU, but is also crucially 
dependent on the assumed magnetic field strength , density , and 
pressure at distances of 200 to 6 OO AU; reexpansion of the compressed 
plasma is enhanced for higher Alfven and sonic speeds. This expan- 
sion results in a positive correlation between density and radial 
velocity in the solar wind If only the region between forward and 
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reverse corotating shocks Is considered. At large radial distances 
this region expands sufficiently so that there is a net positive 
correlation in the solar wind as a whole; a resultant decrease In 
both the average radial velocity and r^n at Increasing distances 
occurs as a result of the increasing correlation. The decrease in 
amplitude of radial velocity fluctuations does not imply large 
Increases In proton temperature as considerable adiabatic compression 
can occur while actually decreasing the average solar-wind tem- 
perature; an estimate based on a crude treatment of the electron 
thermal properties suggests a plateau between 3 to 6 AU with maximum 

average- ion temperatures on the order of 8x10^ K. 

We have also demonstrated the rather surprising conclusion that 
at large radial distances the average solar-wind azimuthal velocity 
may ' be in a direction opposed to the solar rotation. Assuming 
a solar— wind tangential velocity of 8.5 km/sec at 1 AU, values of 


2.5 km/sec at 4 AU and 1.2 km/sec at 6 AU would be expected in the 
absence of stream-stream interactions. In fact, our calculations 
show that the same mechanism that reduces the radial solar-wind velo- 


city also results 


in azimuthal velocities of 0.25 km/sec at ^ AU and 
-1.7 km/sec at $ AU (Fig. 5) in corotatijig models*, for spherical ly- 
symmetric models, a velocity of +0.2 km/sec at 6 AU is expected. 

The total difference at 6 AU due to density-velocity anticorrelation 
is between 1 to 3 km/sec which corresponds to angular deviations of 
from 0.1 to 0.4®. Measuring such relative changes in angle between 
1 and 6 AU is well within the capabilites of present 
plasma experiments; if care is taken in determining 


'*t^ 1 1 "TAifWijMi iT^’iri'rt'ii'rari-iwiir iniM-Wfwr 1 ~ ' 
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detector and spacecraft alignments, an absolute measurement of the 
counterrotation should also be possible. Continuous coverage, 
however, is an Important requirement for such studies, as large 
variations of occur due to stream interaction regions, Alfven 
waves, and other sources of variation. 

We have also predicted that the major fraction of angular- 

momentum flux will be carried by magnetic fluctuations (large scale 

MHD fast waves after Schubert and Coleman, 1968) at distances on the 

order of 3 to ^ AU. This effect should also be observable, despite 

some controversy that has recently been attached to an apparent r“^*^ 

dependence of the azimuthal magnetic field. Parker and Jokipii (1976) 

have recently demonstrated that variations in<V>can account for the 

anomalous radial gradient observed on Pioneer 10; the same correction 

should be made for when computing radial gradients of angular 

momentum flux. An additional problem that impacts estimates of 

radial gradients of field components ts correcting for the reverse 

in direction due to sector structure (Rosenberg £l_. , 1975); as 

angular-momentum flux depends only upon the product of B and B. 

r (|) 

this problem will not influence observations of momentum flux. 
Determinations of radial gradients in angular-momentum flux, like 
other solar-wind properties, can best be made by comparison between 
data obtained simultaneously from two spacecraft at different radial 
distances; otherwise temporal variations (e.g. , solar-cycle velocity 
variations) might bias the observat ions . 
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Average magnetic field components are also affected by stream 
interactions; if a correlation between and at the Inner 
boundary is assumed, a 10% decrease in <rB^> between 1 and 6 AU 
occurs. However, such a correlation is not especially large In the 
solar wind (Goldstein and Siscoe, 1972); several other models tested 
without an initial correlation showed decrease in<rB^>of between 
0 and 10%. 

In summary, nonlinear fluctuations are a significant effect in 
determining the radial gradients of the solar wind at distances as 
small as 0.2 AU; at distances greater than 1 AU nonlinear fluctua- 
tions dominate the behavior of the radial gradients. Therefore, it 
would be desirable that future studies of solar wind gradients 
should directly average .conserved quantities (left hand side of 
eqs. 6- 12), as well as the usual parameters such as magnetic field, 
density, velocity, temperature, and so on. Such averages will allow 
direct estimates of the importance of nonlinear fluctuations in 
determining radial gradients, and are essential for investigating 
effects that are ignored by conventional MHO models, such as energy 
or momentum loss to nonthermal particles or newly ionized solar wind 
neutrals. 
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figure CAPTIONS 

Figure 1. The radial dependence of the average value of with an 
initial densi ty-velocl ty anticorrelation (Case 1, thick 
lines) and without an initial correlation (Case 2, thin 
lines). Spherical ly-symmetric calculations for Cases 1 
and 2 are shown as solid lines, corotating calculations 
as dashed lines. A spherical ly~symmetrlc calculation for 
a large density-velocity anticorrelation (Case 3) is shown as 
alternating long and short dashes. 

Figure 2. The radial dependence of solar-wind number density and 

the product of density and velocity. The long-short dashed line 
is, the product pf average density multiplied by' and the average 
product of average density multiplied by r^ and the average 
velocity. The short dashed line is the difference between 
the number density flux, a conserved quantity, and the 
product just described. 

Figure 3. Terms entering the radial momentum flux equation, all 
quantities have been divided by the proton mass and 
multipl led by r^. 

Figure The ion temperature as a function of radial distance. The short- 
long dashed curve is the electron temperature that results 
from adiabatic changes (y = 5/3)* The (solid, dotted, 
dashed) line represents the ion temperature that results if 
a polytrope law for the electron gas is used with index 
of (1.2, 1.5, 5/3). 
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Figure 5. The radial dependence of the average value of rV^ with 
an initial density-velocity anticorrelation (Case 1, 
thick lines) and without an initial correlation (Case 2, 
thin lines). Spherical ly-synmetric calculations are 
shown as solid lines and corotating calculations as 
dashed lines. 

Figure 6. Angular-momentum flux for the spherically symmetric model 

with an initial density-velocity anticorrelation (Case 1). 

The heavy solid line is the quantity r^<nV^V> ; the 

heavy dashed line is the quantity r^(<nV^V^> - <n><V^xV^>) . 

The thin solid line is the quantity -r^<B^B^>/ypmp; the 

thin dashed line is -r^(<B B > - <B ><B >)/y m . Units 

r (}> r (j> Op 

are AU^ (protons/m^) (m/sec) 

Figure 7. Angular- moment urn flux for the corotating model; variables 
shown as for Figure 6. 

Figure 8. The radial dependence of the quantity " 

<V 'xB >) = A, for a corotating calculation with an Ini- 
4> r 1 

tial correlation between magnetic field and radial 
velocity assumed (Case 4). A 1/r power law is shown as a 
solid line, while the best fit power law dependence of 
1/rl. 05 between 1 and 6 AU is shown as a dashed line. 

The dots are points from the numerical solution. Units 
are gamma X m/sec. 
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